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Abstract. We study the validity of cosmic distance duality relation between angular diameter 
and luminosity distances. To test this duality relation we use the latest Union2 Supernovae 
Type la (SNe la) data for estimating the luminosity distance. The estimation of angular 
diameter distance comes from the samples of galaxy clusters (real and mock) and FRIIb 
radio galaxies. We parameterize the distance duality relation as a function of redshift in six 
different ways. Our results rule out some of the parameterizations significantly. 
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1 Introduction 

To understand the physics behind the late time cosmic acceleration is one of the most chal- 
lenging problem in the cosmology. Many observational methods using distant type la super- 
novae (SNe la), Cosmic Microwave Background (CMB), Galaxy Clusters, Baryon Acoustic 
Oscillations (BAO) etc. have been used to probe the mechanism behind this positive ac- 
celeration [1]. These tools are based on the fundamental assumption that total number of 
photons are conserved on cosmic scales. This is one of the key assumptions in the relationship 
between the luminosity distance, d,L, and the angular diameter distance, d,A, [2], i.e. 

dL 1 • (1-1) 



d A (i + z y 

This equation is known as Distance Duality (DD) relation [3]. It is natural to check the 
validity of this equation which may play a key role in observational cosmology. In particular 
it plays a vital role in galaxy observations, CMB observations and gravitational lensing 
[4]. Furthermore, violation in the DD relation may point to a failure of the metric theory of 
gravity in explaining the background dynamics of the universe and emergence of new physics. 
Whether this DD relation can shed some light on the presence of exotic physics or not was 
first explored by Bassett & Kunz (2004) [2]. They ruled out non-accelerating models of 
universe (replenishing dust model) by more than 4 — a level. In this context to study the 
validity of DD relation (1.1), we analyse the following red-shift dependence of DD 

Uzan et al.(2004) used the combined measurements of Sunyaev-Zeldovich effect and 
X-ray emission data of galaxy clusters to study the violation in the DD relation [5]. They 
showed that if this relation does not hold then the angular diameter distance measured from 
the clusters is d c j^ steT {z) = cIat] 2 , and hence the DD relation (1.2) gets modified to 

r)(z) = — . (1.3) 

They found no significant violation of DD relation as they obtained the value of r\ close to 
1. Their data set consisted of 18 galaxy clusters. Further, De Bernardis et al.(2006) found 
no departure from DD relation i.e. r] = 1 at la C. L.[6]. They used bigger sample of 38 
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galaxy clusters in the redshift range of 0.14 < z < 0.89. Later on Corasaniti (2006) modeled 
the intergalactic dust in terms of the star formation history of the universe and forecasted a 
deviation in the DD relation due to the presence of the cosmic dust extinction [7]. 

Recently Holanda et al.(2011) used two different data sets of galaxy clusters to check the 
consistency of DD relation [8] . These two data sets were based on two different geometrical 
descriptions of clusters. By assuming two parameterizations for rj(z), they showed that the 
data set which is based on the elliptical /3 model of clusters gives better fit with DD relation 
as compared to the data set based on spherical (3 model. Furthermore, Holanda et al. (2010) 
checked the validity of DD relation using the galaxy cluster sample and SNe la with different 
parameterizations of rj(z) [9]. Hence they obtained bounds on the parameters of rj(z). They 
concluded that the best fit values of the parameters of the 77(21) parametrization obtained 
through the data set based on spherical f3 model of cluster are not consistent with the DD 
relation. 

Lazkoz et al. found no evidence of violation in this relation at 2a level [10]. They used 
SNe la and CMB + BAO as standard candles and standard rulers respectively to test the 
validity of this relation. 

Further Avgoustidis et al. demonstrated that the DD relationship can be used to put 
constraints on the cosmic transparency [11]. They obtained the distances by using SNe la 
data, BAO and measurements of the Hubble parameter, H(z). Following the same line of 
thought More et al. also explored the cosmic transparency (conservation of photons phase- 
space density) by using distance measures [12]. 

More recently Lampeitl et al. analyzed the DD relation and found no proof of violation 
at the one sigma level. For this analysis, they used SDSS II SNe data and local BAO 
measurements at redshift z= 0.2 and 0.35 [13]. 

The aim of this paper is to reanalyze the validity of DD relation in a more comprehensive 
manner by using different data samples and parameterizations. In this work we have used 
three different data-sets of galaxy clusters, radio galaxies and a Mock data set to determine 
angular diameter distance, along with Union2 sample of SNe la. We have analyzed the 
validity of DD relation by using six different parameterizations of tj(z). 

The plan of the paper is as follows. In the next section we describe various parame- 
terizations considered in this work. The details of the data and methodology used is also 
explained in this section. We conclude with a section on results and discussion. 

2 Parameterizations, data and method 
2.1 ??(z) parameterizations 

In this work we parameterized the rj(z) with one index and two index parameterizations. 
These relations are inspired by model independent parameterizations for the dark energy 
equation of state [14]. We are parameterizing 77(21) whose value stays one when photon 
number is conserved, gravity is described by a metric theory and photons travel on null 
geodesies. Any significant violation from the DD relation will hint at the emergence of new 
physics. Since all these assumptions are reasonably well tested [2], we expect 77 to stay 
close to unity. To model any departure from unity we parameterize 77 with six parametric 
representations for a possible redshift dependence of the distance duality relation: 
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Two index parameterizations 



T]l(z) = 770 + T)\Z 

mi( z ) = m + m - 



l + z 



z 



(2.1) 
(2.2) 



rjni(z) = m + ? ?5 



z 



(2.3) 



mv{z) =T}Q-r) 7 



(l + z) 2 
ln(l + z) 



(2.4) 



• One index parameterizations 



rjvi(z) = 779 




exp 



( 



l + z 



z 



) 



(2.6) 



(2.5) 



Both, one and two index parameterizations have their limitations and advantages. The 
two index parameterizations are more flexible in comparison to one index parameterizations 
which may completely dominate over data. 



We have used three different data sets to calculate the angular diameter distance. For the 
luminosity distances, we use latest Union2 SNe la data [15]. Now we consider that pair of 
galaxy cluster /radio galaxy and SNe for which Az < 0.005 [9]. Because of this condition, the 
number of data points are limited to 24, 222 and 12 for the data set I, II and III respectively. 
If there are more than one supernova satisfying Az < 0.005 for a given cluster /radio galaxy 
we choose the nearest SNe to the given cluster /radio galaxy. In case of multiple SNe at same 
redshift value for a given cluster /radio galaxy we choose that SNe which has less error bars. 

• Data Set I: This sample consists of 25 galaxy clusters (isothermal elliptical f5 model, 
with concordance model for the cosmological distance-redshift relationship) [16]. By 
combining together the Sunyaev-Zeldovich temperature decrements and X-ray surface 
brightness observations, one can obtain the angular diameter distance for clusters. The 
redshift interval for clusters in this sample is 0.02 < z < 0.78. After applying the 
selection criterion, the number of data points reduce to 24 in this sample. 

• Data Set II: This bigger data set contains 578 mock values of angular diameter distances 
of mock clusters, and for 222 mock clusters we have SNe la satisfying Az < 0.005. 
This catalog is created by assuming fiducial cosmology from the WMAP7 + BAO + 
Hq results [17]. This data set assumes spherical isothermal j3 model for Intra-cluster 
medium (ICM). This mock catalog is created for overlap of ACT/SPT + eROSITA with 
25% Gaussian scatter error which are similar to the recent measurements of angular 
diameter distance [18]. The estimate of angular diameter distances strongly depends 
on the ICM model. The redshift distribution of this sample is from 0.05 < z < 0.76. 
For more details see the ref.'s [19, 20]. 



2.2 Data 
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• Data Set III: In this sample the calculation of angular diameter distance is obtained by 
using the physical size of extended radio galaxies [21]. This data set contains 20 radio 
galaxies up to redshift z = 1.8 . For 12 galaxies we could find corresponding SNe la 
where the condition Az < 0.005 is satisfied. 



2.3 Method 

We perform the x 2 analysis to fit the parameters of the assumed parameterizations. 

Where rjth is the assumed form of parameterizations given by eq.(2.1) to eq.(2.6) and rj b s 
is the observed value of 77, which is calculated by using di and (1a at a particular value of 
redshift. The unknown parameters in the parameterizations are denoted by the variable p. 
To draw the likelihood contours at 1, 2 and 3a, Ax 2 = X 2 ~ Xmin = 2 -3, 6.17 and 11.8 
respectively in two dimensional parametric space. For one parameter fitting, A% 2 = 1,4 and 
9 for 1, 2 and 3 a C.L. respectively. 

In order to find rj b s for cluster, we need the observed value of d c j^ ster and cLl as men- 
tioned in eq. (1.3). The luminosity distance is obtained by using the latest Union2 SNe la 
data. The angular diameter distance, d c j^ ster , is obtained directly by using the data set I 
and II. While calculating the error bars on r] Q b s , we only considered the errors in the angular 
diameter distance values. The errors in the di values are quite small compared to the errors 
in dA, and hence are neglected [9]. We know that the variance of a dependent variable f(xj), 
that depends on, say, n number of independent variables Xj, is given by the error propagation 
equation as 

Of 



dxj 



(2. 



Hence, the error in rj b s for galaxy clusters (mock and real) data set is 



(1 + Zjf 

Similarly the rj b s for the radio galaxies can be obtained by assuming d A adw = dA (R» A. 
Daly, private communication) in eq. (1.2). The angular diameter distance for radio galaxies 
is given in terms of the dimensionless coordinate distance, y, as 

where y = aorHo, and a$r is the coordinate distance [21]. Treating Hubble constant, Ho, 
as a nuisance parameter we marginalize over it by assuming gaussian prior, Hq = 74.2 ± 
3.6 kms~ l Mpc~ l [22]. Finally, the uncertainty on the calculated value of rj b s , for radio 
galaxies can be obtained by using the formula 

<W) - a lA(i) [dA(i) {1 + z . )]2 ■ I 2 ' 11 ) 

In this analysis distances are obtained by assuming the flat A CDM universe. 
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3 Results and Discussions 

In this work, we study the validity of DD relation using different data sets. We assume six 
general parameterizations of rj(z) completely in analogy with the varying equation of state of 
dark energy, cu(z). In these different parameterizations it is assumed that r](z) evolves with 
redshift. 

Our results are summarized as follows: 

1. We find the best fit values of the parameters in all the six parameterizations with 
data set I. This data set contains angular diameter of 24 galaxy clusters. Using \ 2 
minimization technique we obtain the best fit values of the parameters within la range 
shown in Table 1. 



Table 1. Best fit values for all parameterizations - data set I 







Parameters 




Parameters 


1.217 




= 0.999 ± 0.144 


m 


= -0.058 ±0.507 


1.216 


m 


= 1.007 ±0.170 


m 


= -0.118 ±0.822 


1.216 


m 


= 1.013 ±0.205 


Vs 


= -0.198 ± 1.345 


1.216 




= 1.003 ±0.157 


m 


= 0.087 ±0.651 


1.337 


m 


= 1.194 ±0.061 






1.169 




= 1.01 ±0.051 







Here xt is reduced-^ 2 , or x 2 P er degree of freedom. All the two index parameterizations 
do support DD relation with in 1 a C.L. as shown in Fig. 1 for data set I. In one index 
parameterizations, ijy(z), shows significant deviation from DD relation, where as rjvi(z) 
is in agreement with the DD relation (see Fig. 2). 

2. The data set II which contains the mock catalog of angular diameter distances of 
galaxy clusters is the biggest data used for this purpose so far. It consists of 222 galaxy 
clusters after using the selection criterion (Az < 0.005). The best fit values for all the 
parameterizations with in la are shown in Table 2. 



Table 2. Best fit values for all parameterizations - data set II 



Xt/ 




Parameters 




Parameters 


1.076 


Vo 


= 0.973 ± 0.048 


m 


= -0.108 ±0.159 


1.078 


m 


= 0.971 ± 0.058 


m 


= -0.138 ±0.267 


1.080 


m 


= 0.972 ± 0.072 


m 


= -0.141 ±0.453 


1.077 




= 0.972 ± 0.053 


V7 


= 0.124 ±0.208 


1.228 


m 


= 1.163 ±0.019 






1.070 


V9 


= 0.972 ±0.016 







As shown in Fig. 3, all two index parametrization show deviation from DD relation 
at 3a level. While Fig. 4, which shows the likelihood plot for r]y and rjyj, again 
indicates that in case of r\y there is a significant deviation from DD relation. Similarly 
rjyj only marginally accommodates the distance duality relation. Thus, the one index 
parameterizations also do not support the DD relation convincingly. It is important 
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Figure 1. la, 2a and 3cr contours in r]i — r)j plane with data set I. The position of filled circle in the 
contours indicate the point where r\(z = 0) = f . The position of empty circle indicates the best fit 
value of the parameters. 




Figure 2. Likelihood distribution function from the data Set I. The position of empty circle indicates 
the best fit value of the parameters 



to note that five out of six parameterizations do not give a substantial support to 
the DD relation with this bigger mock data set. It may indicate that the underlying 
assumptions used to create this mock data set need revision. 

3. In the data set III, the angular diameter distance is obtained by using the extended 
physical size of radio galaxies. After using the selection criterion (Az < 0.005) we are 
left with 12 data points in this set. The best fit values of the parameters using % 2 
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Figure 3. la, 2a and 3a contours in m — r)j plane with data set II. The position of filled circle in the 
contours indicate the point where r\(z = 0) = 1. The position of empty circle indicates the best fit 
value of the parameters. 




Figure 4. Likelihood distribution from Mock data sample of galaxy clusters ( data Set II ). The 
position of empty circle indicate the best fit value of the parameters. 



minimum method at la level are given in Table 3. 

As shown in Fig. 5 and 6, this data set does not support DD relation for both 2 index 
and 1 index parameterizations convincingly. 

4. It is important to note that all two index parameterizations show degenerate behaviour 
with the given data sets. These parameterizations are in complete concordance with 
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Table 3. Best fit values for all parameterizations - data set III 



Xz/ 




Parameters 




Parameters 


0.944 


Vo 


= 1.063 ±0.198 


Vi 


= -0.180 ±0.244 


0.971 


V2 


= 1.099 ±0.274 




= -0.415 ±0.632 


1.021 




= 1.099 ±0.385 


m 


= -0.749 ± 1.624 


0.958 




= 1.081 ±0.235 




= 0.282 ± 0.404 


1.592 


m 


= 1.564 ±0.070 






0.870 


V9 


= 1.059 ±0.047 








Figure 5. la, 2 a and 3er contours in rji — rjj plane with data set III. The position of filled circle in 
the contours indicate the point where r\(z = 0) = 1. The position of empty circle indicates the best 
fit value of the parameters. 



DD relation within 1 a level for data set I, shows deviation at 3 a level for data set II, 
and in agreement with DD relation within 2 a level for data set III. So the behaviour 
of all two index r](z)'s strongly depends upon the data sets chosen here. 

5. The Xu f° r W parametrization with data sets I, II and III is approximately 1.34, 1.23 
and 1.59, respectively. Hence the one index parametrization, rjy, shows significant devi- 
ation with all the three data sets. The rjvi parametrization marginally accommodates 
the DD relation with data set II and III, and convincingly accommodates data set 
I. It is evident that for one index parameterizations the chosen parametrization prior 
dominates over the data sets. Therefore the one index parameterizations analysed here 
are not a good choice for analysing the DD relation. This is in contrast with the be- 
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Figure 6. Likelihood distribution function from data set III. The position of empty circle indicate 
the best fit value of the parameters. 




Figure 7. Dashed, solid and filled contours correspond to data set I, data set III and combined data 
set (I + III Respectively. 



haviour shown by two index parameterizations in which data seems to dominate. To 
understand behaviour of the parameterizations in a qualitative manner we plot them 
using their best fit values along with real data sets (Figs. 8 and 9). 
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Figure 8. Variation of parameterizations with redshift using best fit values given in Table 1 . The 
points with error bars are from data set I. Solid line correspond to two index parameterizations. Dot 
and dashed lines correspond to one index parameterizations. 




Figure 9. Variation of parameterizations with redshift using best fit values given in Table 3. The 
points with error bars are from data set III. Solid lines correspond to two index parameterizations. 
Dot and dashed lines correspond to one index parameterizations. 



It is clear from the figures 8 and 9 that parametrization r]y{z) is not a good choice 
for modeling the DD relation. It shows a very steep deviation from 1 even at very 
low redshift (z ~ 0), where we expect DD relation to hold. Hence r\y parametrization 
is clearly dominating over the data sets. But on the other side all the two index 
parameterizations stay close to unity (see Figs. 8 and 9). 

6. Out of these three data sets, as expected, the bigger data set of mock galaxy clusters 
(data set II) gives tighter constraints on various parameterizations. But most of the 



parameterizations show significant deviation from DD relation. It is important to note 
that the mock galaxy clusters data is generated by assuming the spherical isothermal 
/3 model for clusters. But the real galaxy cluster data (data set I) which is obtained by 
assuming isothermal elliptical f3 for clusters is in good agreement with DD relation. 

7. The main advantage of using the radio galaxy data is shown in Fig. 7 where the contours 
(dashed lines) correspond to radio galaxies are aligned to contours corresponding to 
clusters (solid line). We also show the contours by using the combined data set of real 
galaxy clusters and radio galaxies. As expected the constraints on the parameters are 
improved. 

8. Overall, the pattern seen between one and two-parameter parameterizations suggest 
that for the two-parameter case the data dominates the result, whereas for the one- 
parameter case the chosen parameterization prior dominates the result. 

When we finished this work, we came across a paper by Li et al. who also have checked 
the violation of DD [23]. However, we have used more data sets, and we have worked with 
six parameterizations. 
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